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ABSTRACT 

We have detected and confirmed five water maser complexes in the Andromeda Galaxy (M31) using the 
Green Bank Telescope. These masers will provide the high brightness temperature point sources needed for 
proper motion studies of M31, enabling measurement of its full three-dimensional velocity vector and its geo- 
metric distance via proper rotation. The motion of M3 1 is the keystone of Local Group dynamics and a gateway 
to the dark matter profiles of galaxies in general. Our survey for water masers selected 206 luminous compact 
24 /im-emitting regions in M31 and was sensitive enough to detect any maser useful for '^10 /ias yr"' astrom- 
etry. The newly discovered masers span the isotropic luminosity range 0.3-1.9 xIQ'^Lq in single spectral 
components and are analogous to luminous Galactic masers. The masers are distributed around the molecular 
ring, including locations close to the major and minor axes, which is nearly ideal for proper motion studies. We 
find no correlation between 24 ^m luminosity and water maser luminosity, suggesting that while water masers 
arise in star-forming regions, the nonlinear amplification pathways and beamed nature of the water masers 
means that they are not predictable based on IR luminosity alone. This suggests that there are additional bright 
masers to be found in M3 1 . We predict that the geometric distance and systemic proper motion of M3 1 can 
be measured in 2-3 years with current facilities. A "moving cluster" observation of diverging masers as M31 
approaches the Galaxy may be possible in the long term. 

Subject headings: astrometry — galaxies: individual (M31) — galaxies: ISM — Local Group — masers — 
radio lines: galaxies 



1. INTRODUCTION 

The dominant galaxies in the Local Group, the Milky Way 
and Andromeda (M3 1 ), are likely to collide in the next 5-10 
Gyr (iLoeb et al.ll2005h . In the meantime, however, we would 
like to know the proper (transverse) motion of M3 1 with re- 
spect to the Milky Way. Peebles et al. (2001) show that mea- 
surements of the transverse motions of Local Group galaxies 
will critically test the fundamental assumption that luminosi- 
ties trace the mass of galaxies. The proper motion of M31 
is key to our understanding of the Local Group's future (and 
past) dynamical evolution as we ll as the density profiles and 
distrib ution of dark matter halos (ILoeb et al.ll2005l: iReid et al] 
I20T0I). 

ILoeb et all (120051) suggest that GAIA, SIM, or the Square 
Kilometer Array will provide the ^10-20 ^as yr"' astrom- 
etry r equired for a proper motion determ ination of M31 
(e.g., Ivan der Marel & Guhathakurtal 120081) . But Galactic 
analog masers in M31 observed with sensitive Very Long 
Baseline Interferometers (VLBI) over a time baseline of 
a few years would be adequate. Indeed, 22 GHz wa- 
ter (H2O) masers h ave been discove r ed in M31's satel- 
lite galaxies M33 (IChurch well et all 119771) and IC 10 
(iHenkel Wouterloot. & Ballvl 119861) 7^ and iBrunthaler et al.l 
(I2OO5L I2007r have measured their proper motions, yielding 
their 3-dimensiqnal vel ocities with respect to the Milky Way. 
IBrunthaler et al.l (120051) also measured the proper rotation of 
M33, which provides a geometric determination of its dis- 
tance via "rotational parallax." Loeb et al. (2005) show that 
a backwards-time-evolution model has M33 colliding with 
M31 in the past. Since M33 does not bear the morphological 
signatures of such an interaction (but does show star forma- 
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tion possibly triggered by a close pass), one must assume that 
M3 1 was elsewhere at the time and thus has a non-negligible 
tangential motion, of order 100 km s~' dLoeb et al. 2005|; 
IBrunthaler et al.ll2007[ but see Ivan der Marel & Guhathakurtal 
(120081) ). 

Detecting ma s ers in M31 has histori cally been problematic 
(lSullivanlll973l) . iGreenhill et al.l ([T995h and llmai et all (l200Tl) 
conducted pointed H2O maser surveys toward H II regions 
in M31, reaching Icr rms noise levels of 29 and 70 mJy, re- 
spectively. Claussen & Beasley (private communication) con- 
ducted a survey for H2O masers in the nuclear region and 
across most of the molecular ring with the Very Large Array, 
makin g no detections and reachi ng an rms noise of 30 mJy per 
beam. ISjouwerman et all (120101) have detected the first maser 
of any type in M31, a Class II methanol (CH3OH) maser at 
6.7 GHz, but no associated water masers were detected. The 
water maser non-detections to date are all consistent with the 
aggregate star formation rate of M31, and these surveys were 
not sufficiently sensitive to have detected a typical Galactic 
maser at the distance of M31 (Brunthaler et al. 2006). The 
lack of Galactic analog water maser detections in M31 has 
simply been a problem of sensitivity and of locating likely 
sites of maser activity, both of which now have remedies. 

Galactic analog H2O masers associated with star formation 
are detectable in M3 1 in reasonable integration time using the 
Green Bank Telescop43 (GBT); they akeady have been de- 
tected in significantly more distant galaxies (Darling, Brogan, 
& Johnson 2008). But a full survey of the star-forming re- 
gions of M31 with the sensitivity to detect Galactic analog 
water masers is prohibitive because the galaxy is so large on 
the sky. We thus combined an order of magnitude sensitiv- 
ity improvement compared to the best previous survey with 

^ The National Radio Astronomy Observatory is a facility of the National 
Science Foundation operated under cooperative agreement by Associated 
Universities, Inc. 
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Fig. 1. — Spitzer 24 map of M31 jGordon et alj|2008b . The circles indicate tlie 206 pointing centers (tlie circles are to scale, showi ng the 33" FWHM 
beam). Crosses mai'k the detected masers and are enlarged for clarity. The orange cross indicates the 6.7 GHz methanol maser detected bv lSiouwerman et al.l 
120101) . Colors indicate the systemic (green), red, and blue spectrometer tuning centers at —300, —100, and -500 km s"' , respectively. 



a physically motivated pointing guide to the most likely sites 
of water maser activity. In this Letter, we describe the source 
selection, the observations and data reduction, and the newly 
discovered water masers. We also discuss the prospects for 
using these masers for proper motion studies of M3 1 . 

Throughout this Letter, we use heliocentric velocities (op- 
tical definition) and assume th e following about M3 1 : a s ys- 
temic velocity of -30 km s~ ' ("de Vaucouleurs et al.lll991h . a 
distance of 780 kpc (McConn achie et al. 2005), a major axis 
position angle of 38° (E of N), a constant inclination of 77°, 
and a central coordinates 00:42:44.3, +41:16:09 (J2000). 

2. SAMPLE SELECTION 



A Spitzer 24 /im map of M31 (iGordon et al.ll2008h guided 
the GBT observations in two important ways: source selec- 
tion and completeness. The association of water masers with 
(ultra)compact H II regions is well known, but the molecular 
gas as traced by CO emission does not necessarily indicate 
these regions. Likewise, Ha-selected H II regions may select 
against dust-enshrouded ultra-compact H II regions. Compact 
24 yum emission, however, is a good indicator of likely H2O 
masers, and the rough relationship between far-I R emission 
and water maser luminosity is well-known (e.g., iJaffe et all 
119811: ICastangia et alj|2008[) . We therefore selected compact 
(unresolved at the 23.5 pc resolution of the Spitzer map) 
24 /im sources associated with the dusty molecular — and 
presumably star-forming — regions in M3 1 . We constructed 
a catalog of 206 of these objects working from the brightest 
down to the point where most of the 24 fim. emission becomes 
extended, about 4 MJy steradian"' at peak specific intensity. 
Although this pointed survey does not include all star forma- 
tion in M31, it does include a large fraction of the likely sites 
for strong H2O maser activity, and by working down in lumi- 
nosity through the catalog of compact 24 /im sources, we have 
surveyed a large fraction of the total ongoing star formation 
in the galaxy (e.g., Tabatabaei & Berkhuijsen 2010) and thus 
most of the likely H20-bearing regions (^. Figure [T] shows 
the 24 /im map of M3 1 with the pointing centers and primary 
beam size for our GBT water maser survey. 

3. OBSERVATIONS AND DATA REDUCTION 

We observed the 6i6-523 22.23508 GHz ortho-water maser 
line toward 206 24 /^m sources in M31 with the GBT in 2010 



October through December A 670 km s ' (50 MHz) bandpass 
was centered on a heliocentric velocity of -300 km s"' for 97 
sources in the central parts of the galaxy and along the minor 
axis, on - 100 km s"' for 52 sources in the redshifted northeast 
wedge of the galaxy, and on -500 km s"' for 57 sources in the 
blueshifted southwest wedge of the galaxy (Table[Tli. 

Observations were conducted with the dual K-band re- 
ceivers in a nodding mode in two circular polarizations with 
0.16 km s~' (12.2 kHz) channels and 9-level sampling. The 
time on-source was 5 minutes except for sources that were re- 
observed to confirm or refute possible lines (typically 10 min- 
utes; see Table [T]i- A winking calibration diode and hourly 
atmospheric opacity estimates were used for flux density cal- 
ibration. Opacities ranged from 0.05 to 0.13 nepers but were 
typically 0.07 nepers. The estimated uncertainty in the flux 
density calibration is ^20%. Pointing was typically good to 
within a few arcseconds and the largest pointing drifts dur- 
ing observations were no m ore than 6" . The res olution of the 
24 /im Spitzer image is 6" ( iGordon et al.ll2008l) . so the unre- 
solved IR sources remained within the 33" GBT beam even 
during the largest pointing drifts. The 33" beam (FWHM) at 
22 GHz spans 125 pc in M31. 

After averaging polarizations, spectra were Hanning 
smoothed and subsequently Gaussian smoothed to obtain a 
final spectral resolution of 3.3 km s~' (244 kHz) for the line 
search. Polynomial baselines, typically of fifth order, were 
fit and subtracted to obtain flat and generally uniform-noise 
spectra. Spectral rms noise measurements for each pointing 
center are listed in Table [T] Non-detection spectra generally 
did not show any features greater than 3ct, corresponding to 
roughly 10 mJy, at lea st an order of magnitude more sensitive 
than previous surveys (iGreenhill et al.lll995tllrnaretani200Tl 
Claussen & Beasley, priv. comm.) and likely the weakest line 
that can currently be used for VLBI proper motion studies. 
All data reduction was performed in GBTIDLQ. 

4. RESULTS 

We have detected water maser emission in five out of 206 
24 /im-selected regions in M31 (Table [T]i. Most of these are 

' GBTIDL jhttp : / /gbt idl ■ nrao ■ edu/) is the data reduction pack- 
age produced by NRAO and written in the IDL language for the reduction of 
GBT data. 



Water Masers in the Andromeda Galaxy 



3 



TABLE 1 
Water Maser Survey Sample 



Source 


Vob. 


rms 


Source 


Vobs 


iTns 


Source 


Vobs 


rms 


Source 




rms 


J2000 


km s"' 


mJy 


J2000 


km s"' 


mJy 


J2000 


km s"' 


mJy 


J2000 


km s"' 


mJy 


003849.2+40255 1 .7 


—500 


3.2 


003906.7+403704.5 


-500 


3.5 


AAOA1A i A f\^^^ C z' 

0039 1 0.2+403725 .6 


-500 


3.2 


A A OA1/' 1 1 /lAO/^IA C 

003916.1+403629.5 


-500 


3.0 


003918.9+402158.4 


-500 


1.8 


003930.2+402 1 06.4 


-500 


3.7 


AATAOO ^ ■ ^A'^'^ 1 C 

003933 .2+4022 15.6 


-500 


3.2 


AA1 A H C 1 /1A1A1 1 C 

003937.5+40201 1.5 


-500 


3.4 




—500 


3.1 


003944.5+402030.4 


—500 


3.3 


AA/lAA/1 ~7 1 A nCO A r\ A 

004004. /+405840.9 


—500 


3.2 


AA/IAIA Ai /lA/inAA 

00403 . 9+4042 30 . 


—500 


3.0 


Uu4U J 1 . 2+40395 2 . U 


—500 


2.9 


r\f\Ar\1 1 "7 ■ Af\A 1 T7 A 

004031.7+404127.0 


—500 


3.3 


AA/IATO /I 1 ylATOC/l 1 

004032.6+403856. 1 


—500 


3.3 


AA/1 ATT 1 1 /I AT f CI 1 

004033.3+403352. 1 


—500 


3.1 


(JU4Uji.o+4Ui24o.O 


—500 


3. 1 


004034. /+403541.2 


—500 


2.8 


004035. 1+403 /Oi. 1 


—500 


4.0 


004035.8+403 /24.6 


—500 


3.6 


/l/l Af\0£i 7i A 1 f\ 1 I 1 ^ 

0040 JO. 1+4101 1 /.J 


—500 


2.1 


00403&. 0+403514.9 


—500 


2.9 


A A /I AlOOi/IA'I/l'JI A 

004038. 6+40343 1 .0 


—500 


3.4 


AA/IAIA A 1 /1AT~7TA C 

004039.4+403/30.5 


—500 


3.4 


004041 .0+40 J IOj.O 


—500 


2.0 


r\f\Af\A'^ 1 1 Af\1ACA C 

004042.1+403454.5 


—500 


3.4 


AA/IA/I"? ") 1 ylA/lTll A 

004043 .3+40432 1 .9 


—500 


3.1 


AA/1 A /IT C 1 /I AT C T A C 

004043.6+403530.5 


—500 


3.3 


(JU4U44.2+4U4440.4 


—500 


3.0 


U0404J. / +40j1j4.j 


—500 


1 n 
1.9 


004046.4+40554 1 .9 


—500 


3.0 


004046 . 5+405 606. 4 


—500 


2.8 


UU405 1 .6+4 1 OUUo . 5 


—500 


2.9 


r\f\A r\c 1 r\ 1 A ac* 1 

00405 1 .9+403249. / 


—500 


3.4 


A A /I AC! Ai /lAniO A 

00405 3 . 0+4032 18.0 


—500 


3.8 


AA/1ACC 1 1 /1AT~7AT 1 

004055 . 1 +403 /03 . 2 


—500 


3.8 


0040 J O.O+404DJO.0 


—500 


2.1 


004056.6+410332.3 


—500 


3.0 


AA/IACA O 1 ylAT/lCI A 

004059.8+403652.4 


—500 


4.8 


A A A 1 A A £1 1 /11ATT/1 A 

004 1 00.6+4 1 0334.0 


—500 


3.7 


004101. o+41040j.o 


—500 


2.9 


0041 02.0+4102 j4.y 


—500 


1.7 


A A /I 1 A"7 i /I AH O 1 1 C 

004 10/. 6+4048 12.5 


—500 


3.5 


AA At 1A A 1 /lA/lA/lA C 

0041 10.4+404949.5 


—500 


3.1 


0041 12.5+410609.7 


—300 


3.4 


004113. 7+403918.6 


-500 


2.2 


0041 13.9+410736.1 


-300 


2.2 


AA A \ t A O 1 At AA1 T 1 

0041 14.8+410923.7 


—300 


3.8 


f\f\ A\\cc^lA^iA^i^\ c 

UU41 15.y+4U4L)l 1 .6 


—500 


3.1 


r\r\A iin 1 • a r\A a 

0041 19.1+40485 /.4 


—500 


3.6 


A A /I 11ACiyll1A/10 

0041 19.5+41 1948.8 


—300 


3.9 


A A /IIIAAi/IIAOII C 

004 1 20.0+4 1 082 1 .5 


—300 


3.0 


UU4 1 2 1 .7 +404947 .7 


—500 


2.4 


004 1 23 .2+405000.6 


—500 


3.4 


A A /I 11/1 Oi/1111C/l ^ 

004124.8+41 1 154.6 


—300 


3.4 


AA .1 1 O C A 1 /lA/lOAA A 

0041 25 .4+404200.4 


—500 


3.7 


004 126.1+404959.1 


—500 


3.7 


f\f\ A -\ '^f c ■ A i^r\/' r\ 

004126.5+41 1206.9 


-300 


3.0 


AA/iiio 1 ■ At im ^ 

004128.1+41 1222.6 


-300 


3.6 


AA ^ 7 TA O 1 Af\Cf\Cf\ C 

004129.8+405059.5 


-500 


2.2 


UU412y. 0+4 1221 1.1 


—300 


3.1 


004130.3+410501. / 


—500 


3.1 


AA/1 1T"7Ai/1AC1/n C 

00413 /. 0+405142. 5 


—500 


3.3 


AA /11TO ^ 1 /1A/1/1A1 O 

004 1 3 8 .6+40440 1 . 2 


—500 


3.4 


004146.7+41 1b46. 6 


—300 


3.3 


f\f\ A -\ AO . ^11 C\f\^ 

004148.2+41 1903.8 


-300 


3.6 


004 1 5 1 . 9+4 1 2442. 1 


-300 


3.3 


AA.,^7£"A A 1 //AC71A O 

0041 59.4+405720.8 


-500 


2.1 


UU42U3. 9+4U49U/. I 


—300 


3.7 


004208. 8+4 1 2639.9 


-300 


3.2 


AA/11AA ^ i T 

004209.5+412832.3 


-300 


3.3 


AA /1111 , A t t AAA A 

00421 1.6+41 1909.4 


-300 


2.7 


004212.3+412415.7 


—300 


2.9 


004213.8+4051 17.7 


—300 


3.9 


004214.8+412508.9 


-300 


2.8 


AA /1110^i/linCl o 

0042 18.7+412751.8 


-300 


2.7 


UU422 1 . /+4 1 2 52 / . 6 


—300 


2.8 


004226.1+410548.2 


—500 


3.8 


004226.4+41281 1.2 


—300 


3.4 


AA-lnA 1 1 /I 1 TAA/I A 

004230. 1+41 2904.0 


—300 


2.8 


004230.3+412935.9 


—300 


4.0 


004230.9+405714.6 


-300 


3.0 


004234.2+413007.3 


-300 


3.0 


AA /11TC /^i /11T1/1AA 

004235.6+413149.0 


-300 


3.2 


004238.6+413150.5 


-300 


3.7 


004238.9+413135.6 


-300 


2.9 


A/'^ /I 1 /I A 1 i At AI H 

004240.1+410222.7 


-300 


3.5 


AA A^ A t Ti/llH/1/^T 

00424 1.3+412246.7 


-300 


4.1 


nr\ A'^ A 1 r\i/ir\cicci 

UU424 1 .y +4U5 155.2 


—300 


3.6 


A A /1 1 /1 1 C i At AAA 1 A 

004242 .5+41 000 1 .4 


—300 


3.0 


004242.5+413155.2 


—300 


3.9 


AA A'~\ A'~\ Ai /11T1CA O 

004242.9+4 13159.8 


—300 


3.7 


(JU4244. 1+413259.2 


—300 


3.4 


A A /I A A A • A 1 1 ^r\0 C 

004244.4+41 1608.5 


—300 


3.7 


004246.2+4101 1 1.4 


—300 


2.4 


AA A'^ AI Ai /11TTTT A 

004247.0+413333.0 


—300 


4.0 


r\r\ A^ A^ c\ , A -t ^ A f\r\ c 

004247.9+413400.5 


-300 


3.1 


AA A^ Ar\ 1 ■ A t -t C\ A C A 

004249 . 1 +4 1 1 945 .9 


-300 


3.1 


004249.3+412507.5 


-300 


3.2 


AA A'^C t Ai/IIT Cm o 

00425 1 .0+4 1 3507.8 


-300 


3.2 


UU4252.J+4I0U14.S 


—300 


2.1 


004252.4+4 1 1 20.7 


-300 


2.8 


AA.71C^ A ■ ^7 ^TTO 7 

004256. 9+413728. 1 


-300 


1.6 


AA A^CO 1 1 A \ f\f\ t C f\ 

004258.2+410015.9 


-300 


3.1 


f\r\ A T r\r\ < A 1 1 ^ c a 

(JU43UU.U+4 13654.2 


—300 


2.5 


A A /I TAI Ci^l'J'71'7 1 

004301.5+413/1 1.2 


-300 


2.8 


r\f\A f AI C i A t 11 A r\ c 

004302.5+413740.5 


—300 


3.5 


AA/1 T r\£L A 1 /I 1 T OAT 1 

004306.9+4 1 3807 . 1 


—300 


3.3 


Uu431 1.3+410459.5 


—300 


3.6 


A A /I 1 1 A 1 i H 1 AAT 1 A 

0043 1 4 . 2+4 1 003 3.9 


—300 


2.4 


r\nA O"^ A liAlAI~\A'7 

004524. 1 +414124. / 


—300 


2.0 


AA A A 3i A 1 A A i O 1 

004d24. J +4 1441 o. / 


—300 


1.9 


UU4324.4+41UoU2.9 


—300 


2.5 


A A /I "! 1 C ^ 1 At AT A^ A 

004325 .6+4 1 0206.4 


—300 


4.4 


AA/1 TT"! I A t t AT~\ 1 

004333.6+41 1432.3 


—300 


2.8 


AA ATI A A 1 At A A C T /ll 

004334.9+4 1 0953 .6 


—300 


3.8 


004 j41./ +414j19.4 


—300 


2.0 


004343.9+4111 J7.0 


-300 


1.8 


AA/1 f/IZ^Ii/ll/l/llO C 

004346.3+414418.5 


—300 


3.8 


AA A'~l A O 1 /1110TA T 

004346.8+41 1239.7 


—300 


2.9 


f\f\ An AO 1 i/llllTTI 

UU4348.1+41 1 133.2 


—300 


3.2 


A A /I 1 An A 1 At 1AC1 O 

004349.4+41 1053.8 


—300 


2.4 


AA/1 1 C A O I A 1 A"! 1 C £. 

004354.8+414/15.6 


—300 


2.9 


AA /ITCCOi/lllOII £. 

004355.8+41 121 1 .6 


—300 


4.1 


004i5o.o+412o29.o 


—100 


2.4 


AA /IOC/" O . A t A o ^ i r 

004356.8+414831.6 


-300 


3.6 


f\f \ A"^ d Ti A \ A O C A A 

004357.7+414854.0 


-300 


3.1 


A A A'~) C O 1 1 At A1^ ^ A 

004358.2+414726.9 


—300 


3.0 


(jU435o. /+414o3 /.5 


—300 


3.0 


(\(\A /I A 1 ^ \ Ay A nnn a 
004401 .5+414909.6 


—300 


3.8 


004403.0+414954. / 


—300 


3.1 


nr\A Af\i c\ 1 A A ^ A ^ A o 
004403.9+413414.8 


—100 


1 A 

3.4 


nr\ A Ar\A n < /iicrti/; i 

Uu44U4.y +4 15016.1 


—300 


2.8 


A A /I A r\/l A ■ /11T7/1C A 

004406.4+4 1 2 /45 .0 


—100 


2.9 


AA/1 A A"7 Ai /linCA T 

00440/. 0+412759. 3 


—100 


3.3 


AA A A CM Oi/IIIIIAT 

00440 / . 8 +4 12110.7 


—300 


3.0 


U044uy.5+411o5o.o 


—300 


1.6 


A A /I 1 A C 1 ^ n An A1 z 

( J044 10.5 +420247 . 5 


—100 


2.9 


AA/1 /111 1 i/llfTlAC 

0044 12.1+4 1 3320.5 


—100 


2.8 


AA/1 A t C 1 1 A t t AAC 1 

004415.3+41 1905.7 


—300 


2.9 


UU4416.U+41495U. / 


—100 


3.1 


0044 1 5 . 2+4 1 3 406 . 6 


—100 


2.8 


AA/1 /I 1 A O i /I 1 1 ATA A 

004419.2+41 1930.9 


—300 


2.9 


A A /1/11AAi/limA1 O 

0044 1 9.9+4 1 220 1 .2 


—300 


3.2 


nr\ AA'^f\^lA\c^f\\ c 

UU442U.2+4 15101.5 


—100 


3.4 


A A /I /1TI Ai HIIACA A 

004423 .0+4 1 2050.9 


—300 


3.3 


AA/1 /in Ti/IITO/H ^ 

004423.3+413842.6 


—100 


2.8 


AA A A">'^ 1 1 A t '~\ ATI T 

004423. /+4l24o 1 .5 


—300 


3.3 


f\f\ A A^ A 

004424.1+4121 17.3 


—300 


2.9 


AA A A ^ A Oiy(1/lA10A 

004424.2+414918.9 


-100 


3.2 


004424.9+413739.1 


-100 


2.1 


f\r\ /i/ii/'ii,iiiAC/i 1 

004426.2+412054.1 


—300 
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Note. — 206 .sources were obsei-ved for 5 minutes except for italicized sources (10 min). Bold sources are HiO masers; 003918+402158, 004343+41 1 137, and 004409+41 1856 
were observed for 15 minutes in two sessions; 004121+404947 was observed for 10 minutes in a single session; and 004430+415154 was observed for 25 minutes in tliree sessions. No 
maser variability was apparent. The rms noise is for 3.3 km s"' channels in all cases. Vbbs is the observed central velocity. 



clearly maser complexes that will likely resolve into individ- 
ual maser spots when imaged interferometrically. Measured 
line properties are listed in Table |2] spectra are plotted in Fig- 
ure |2] The isotropic line luminosities of the individual maser 
components are anal ogous to luminous Galacti c water masers. 
We also confirm the ISjouwerman et al.l (12010b non-detection 
of water at the location of the methanol maser (Figure [T}. 

Wherever there is coverage, the water lines correspond 
to the range of CO (1-0 ) line velocities within the GBT 
beam (iNieten et al\ 120061) . The notable exceptions are 
004409H-411856, in which the second water maser compo- 
nent does not overlap with the observed CO line range, and 
003918H-402158 where there is no CO line coverage. In the 



latter, the observed maser velocities are broa dly consistent 
with the H I rotation curv e and its uncertainties (IChemin et al.l 
l2009l although see ^5^ . 



5. DISCUSSION 

The GBT survey has produced some surprises: (1) the 
strongest water maser complex lies outside of the main star- 
forming molecular ring in M31, and (2) while the 24 /im 
selection method is a success, the 24 /im luminosity is not 
predictive of maser presence or strength. Our survey has re- 
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TABLE 2 
Water Maser Line Properties 



Source 


Line 


(km s"') 


^Peak 

(mJy) 


AVfwhm 
(kms-') 




fSdv 

(mJy km s ' ) 


Lh^o 


V\jjy 1 0+4UZ 1 JO 


1 
i 


—J /4.3U(U.U / J 


J J.j(4. / ) 


1 rifr\ '^\ 
1 V{\J.Z) 


3/(0) 


A AAA^O/ 1 1 A 

U.UUU JZ( 1 i J 




z 


— JDD. /y(U.U4) 


ui. /(4.zj 


1 Q/Tl 1 A 
1 . J(U.i ) 


8/(9) 


A A A 1 O/^l A 




'J 

3 


— JD3.jy(U.US) 


26.y(j.o) 


U. 6(0.2) 


23(8) 


A AAAn/ 1 1 \ 

U.UUU3z(1 1) 




4 


-561.84(0.02) 


122.2(9.1) 


1.13(0.09) 


135(15) 


0.0019(2) 




5 


-559.7(0.4) 


26.0(2.8) 


3.6(1.3) 


99(36) 


0.0014(5) 




6 


-557.63(0.04) 


68.5(7.1) 


0.9(0.1) 


66(11) 


0.00092(15) 




Sum 








447(43) 


0.0063(6) 


004121+404947 


1 


-528.9(0.4) 


17.7(2.4) 


5.1(0.8) 


95(20) 


0.0013(3) 


004343+411137 


1 


-297.15(0.10) 


21.1(4.1) 


1.0(0.2) 


23(7) 


0.00032(10) 




2 


-294.01(0.04) 


44.9(4.8) 


0.74(0.09) 


35(6) 


0.00049(8) 




Sum 








58(9) 


0.00081(13) 


004409+411856 


1 


-245.09(0.07) 


32.7(4.5) 


1.1(0.2) 


38(8) 


0.00053(11) 




2 


-237.4(0.2) 


18.1(3.1) 


2.3(0.5) 


44(12) 


0.00062(17) 




Sum 








82(14) 


0.0011(2) 


004430+415154 


1 


-117.5(0.5) 


9.3(1.4) 


6.8(1.2) 


67(15) 


0.0009(2) 



Note. — Line properties determined from Gaussian fits. Isotropic line luminosities were computed from L/z^o = (0.023 Lq)D~ f Sdv, where D is in Mpc and / Sdv is in Jy km 
IHenkel et aU200^ . Parenthetical values indicate la statistical uncertainties. 
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Fig. 2. — Spectra of the five water maser complexes detected in M31. 
While all of our spectra were smoothed to 3.3 km s"' for the search, the lines 
of 003918+402158, 004343+411137, and 004409+411856 are very strong 
and narrow and we present their spectra with 0.33 km s"' channels. For 
003918+402158, we plot fit residuals offset by -50 mJy. Dotted lines on the 
other spectra indicate the zero point. Numbers below the spectra label the 
fines enumerated in Table |2] The bold horizontal bar s indicate the velocity 
range of CO within the GBT beam iNieten efQf.l2Q06D . 



vealed no correlation between 24 fim luminosity and water 
maser luminosity, suggesting that while water masers arise 
in star-forming regions, the nonlinear amplification pathways 
and beamed nature of the water masers means that they are not 
predictable based on IR luminosity alone. We therefore be- 
lieve that additional strong water masers may be found in M3 1 
because our survey of 24 /im-emitting regions has not been 
exhaustive — it has simply selected the most IR-luminous 
compact sources (Figure[l}. 

It is also certain that not all of the selected 24 /im sources 
are associated compact H II regions: cool stellar atmospheres, 
supernova remnants, planetary nebulae, and background AGN 
and galaxies can all be sources of 24 i^im emission (e.g., 
IVerley et al.l2007HMould et al.ll2008h . Our detection statistics 
among the observed sample are thus not physically meaning- 
ful, although the majority of the 24 /im emitters in the molec- 
ular ring of M3 1 are likely to be compact H II regions. 

5.1. Detection Rate and Star Formation 



The sensitivity of the GBT survey corresponds to a ^200 Jy 
Galactic maser 6 kpc distant. There are 23 H2O masers in 
the Palagi et al. ( 1993) sample and about 30 masers in the 
iBreen et all ( 120101) sample that our survey could detect at the 
distance of M3 1 . Hence, sensitivity is not the only considera- 
tion; the number of water masers detectable in M3 1 will also 
depend on the the substantially lower star formation rate of 
M31 compared to Galactic. 

To estimate the num ber of expected water masers in 
Local Group galaxie s, iBrunthaler et al.l (12006b used the 
iGreenhill et al.l ([r990l) Galactic water maser luminosity func- 
tion (LF) scaled from the Galactic star formation rate to that of 
Local Group galaxies, assuming that maser population is pro- 
portional to the star formation rate. Local Gr oup water maser 
statistics agree fairly well with this approach (IBrunthaler et al.l 
I2OO6I) . We scale the Galactic star formation rate (~4 Mq yr"^ ; 
Diehl et al. 2006) to that of M31 (-0.35 yr"'; Walter- 
bos & Braun 1994, Tabatabaei & Berkhuijsen 2010) to pre- 
dict 2.3 water masers for an isotropic line luminosity limit of 
4 X lO"'* Lq, which is a 3cr line flux limit of 9 mJy and width 
3.3 km s"^. This small number suffers from small number 
statistics, but our five detections suggest that either we have 
detected the majority of bright water masers in M3 1 or that the 
star formation scaling approach is inaccurate. This program 
is the first to detect the tail of the H2O maser LF in M3 1 . 

5.2. Velocities 

A comparison of the maser velocities to the rotation of M3 1 
is of special interest for a geometric distance determination. 
Proper rotation is a measure of angular frequency, and if this 
angular motion corresponds directly to the rotation velocity, 
then a geometric distance may be obtained: D = Vroi / Urm ■ A 
proper rotation distance determination relies on the the masers 
rotating with the galaxy, but they could also exhibit their own 
peculiar proper motions, so a large ensemble of water masers 
would be ideal to reduce this source of error 

In Figure |3] we examine the correspondence between the 
maser velocities and the p rojected and inclined H I rotation 
curve ( IChemin et al.ll2009l) . Since the maser velocities are ex- 
tremely well-determined (sub-km s"'), all of the uncertainty 
in the velocity comparison arises from the innate velocity dis- 
persion of the disk, the statistical uncertainty in the rotation 
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Fig. 3. — Comparison of maser velocities to tlie H I rotation curve of M31 jCliemin et al.||20()9(). Colors match those in Figure[T] except the methanol maser 
is indicated by a green triangle. Left: Observed water maser velocity versus position along the major axis. The dotted lines indicate the expected observed H I 
rotation curves at the radius of each maser. Once the radius is reached, the inclined rotation curve is seen along the major axis. Right: The difference between the 
observed water maser velocities and the projected inclined H I rotation velocity versus position along the major axis. The dashed lines indicate a fiducial rotation 
uncertainty of ±10 km s"' , which is typical of M31 H I rotation curve measurements. The maser velocity uncertainties are much less than 1 km s"' (Table|2). 

curve ('^10 km s"'), and the systematic uncertainties such as 
inclination, wh ich can induce a rota tion curve offset of order 
10-15 km s-i dCorbelli et alJlMol) . While the masers show 
fairly good agreement with the H I rotation, there is a notable 
overall blueshift that may be indicative of preferentially ori- 
ented maser amplification pathways associated with outflows. 
But given the propensity for masers to show peculiar veloci- 
ties, a larger sample is needed to rigorously examine this ef- 
fect. 



5.3. Proper Motions 

The three main contributions to the proper motion of masers 
in M31 are (1) the rotation of the spiral disk, (2) the systemic 
proper motion of M31 itself, and (3) the relative motion of 
masers omitting the galaxy's rotation: 

(1) The direct observation of the rotation of th e M3 1 spiral 
disk is the ( in)famous van Maanen experiment ( Ivan MaanenI 
I1923L Il935h . A rotating disk seen edge-on will show no 
proper motion at the tangent points and the full circular ve- 
locity along the minor axis. M3rs inclined disk has a circu- 
lar velocity at the molecular ring radius of about 250 km s~', 
which becomes a transverse motion of 56 km s"' at the tan- 
gent points or 15 microarcseconds (/ias) year"'. The trans- 
verse motion along the minor axis is expected to be about 
70 ^ as yr~' . Since the rotation curve of M3 1 is well-known 
(e.g.. iChemin etani2009l: ICorbelh et al.llMol) . the measure- 
ment of the proper motions of masers distributed around the 
disk will provide a geometric distance to the galaxy, limited 
roughly equally by the uncertainty in the rotation model and 
in the proper motion measurements. 

(2) The transverse motion of M31 provides the two un- 
known components of its three-dimensional velocity. The 
expected transverse velocity of M3 1 is of order 100 km s"' 
jLoeb et al.. .2005; .Brunthaler et all 120071 but see van der 



Marel & Guhathakurta 2008), or 27 /ias yr"', and measurable 
over two or three years (see below). 

(3) There may be relative motions of masers within M31 
due to peculiar deviations from the bulk rotation of the disk 
and due to the line of sight motion of M31 itself. The for- 
mer peculiar motions are of order 10-20 km s~' (Figure O, 
which is significantly smaller than the velocities discussed 
above and of the same order as the uncertainty in the rota- 
tion curve. These will tend to average out over an ensemble 
of maser complexes. The latter is the moving cluster effect 
caused by M3rs radial velocity toward the Galaxy; M31 gets 
areally larger as it nears, causing relative maser separations to 
grow over time. This effect is small, about 3 fias yr"' of di- 
vergence for the maser pair with the largest separation, 1.79°. 
The size of this effect is equivalent to '^lO km s"', similar 
to maser peculiar motions. It is possible that with long time 
baselines and a rich network of masers one could extract the 
distinct signal of diverging maser spots and thus estimate the 
distance to M3 1 (but with lower precision than the proper ro- 
tation method). 

How long will it take to obtain proper motions, and what 
are the limits on the prope r motions that can be measured? 
The lBrunthaler et all (120051) study of the proper motion of wa- 
ter masers in M33 informs our expectations for M31. Water 
masers appear in many-component complexes that can show 
dramatic variability, including the appearance of new masers 
and the disappearance of others. But over the span of about 
three years, Brunthaler et al. (2005) were able to track several 
persistent maser features and obtain proper motions relative 
to background compact radio sources. Their astrometric un- 
certainties were on average 8-10 iias in each epoch. Proper 
motion uncertainties from a three-year fit were r^5 /ias yr"', 
which would correspond to a transverse velocity of 18 km s"' 
in M31. 
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We estimate that in two 12-hour observations with the High 
Sensitivity Array one can obtain 1 mJy rms noise in 1 km s"' 
channels. This translates into a 6 /xas uncertainty per epoch 
for a 40 mJy line. In a single year we expect a lOcr detection 
of the proper rotation along the minor axis, a 2.5(7 detection 
along the major axis, and a 4.5(7 detection of the systemic 
transverse proper motion of M31 (assuming 100 km s"'). 
Note that proper motion measurement errors improve as t'^l- 
(jReid et al. 2010), so a 2-3 year astrometry campaign should 
be able to determine with high significance both the geometric 
distance to M3 1 and its transverse velocity. 



6. CONCLUSIONS 

We have identified the first H2O masers in the M3 1 based 
on pointed observations of 24 /im-selected regions. The 
24 pm selection method clearly succeeds in identifying water 
masers, but the 24 pm properties are not predictive of maser 
presence or brightness. We suggest that the catalog of water 
masers in M3 1 is incomplete and that an exhaustive survey of 



IR-luminous regions will identify additional masers. 

The newly identified maser complexes are suitable for 
proper motion studies of M31. The next steps for proper 
motion and geometric distance measurements are interfero- 
metric localization and mapping followed by VLBI mapping 
and monitoring. We predict that proper rotation and systemic 
proper motion will be measurable in a few years with current 
facilities, and that the "moving cluster" divergence of maser 
spots as M3 1 approaches the Galaxy may be detectable in the 
long term. 
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